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I. INTRODUCTION 



A. GENERAL 

The present knowledge of rotar) wing unsteady aerodynamics is limited. The 
fundamental closed form solutions by The odorsen [Ref. 1] and Loewy [Ref. 2] provide 
the basis for the theoretical work in this area. While Loewy's work on wake-induced 
flutter helps explain the phenomenon, it also points out the difficulty to be overcome 
which is the closed form theory shows rapid changes in the lift deficiency function 
with changes in reduced frequency (k), wake spacing (h), and frequency ratio ( m ). 

In the past, emphasis in the study of unsteady aerodynamics has tended to focus on 
flutter instability and the effects of unsteady aerodynamics on lift generation and 
torsional loads, and not on performance. In this thesis, the emphasis will be on 
performance and the effect of unsteady aerodynamics on the drag of the airfoil. The 
classic work on this subject was performed by Garrick [Ref. 3]. In this work, Garrick 
found that the drag force could be reversed, acting to propel the airfoil forward. This 
drag reversal is generally referred to as either negative drag or propulsive force. 

The primary objective of this thesis is twofold. First, the effect of rotary wing 
unsteady aerodynamics on the drag of an airfoil will be analyzed. It will be shown 
that this effect can be either substantially greater than that shown by Garrick, or 
substantially less depending on the frequency ratio (w). This change in the effect is 
due to layers of shed vorticity (or wakes) below the rotor blades. Second, with the 
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recent advent of computer panel code- c: >able of calculating the time histories of the 
wakes generated by oscillating airfoil v a dosed form theory is developed for the 
rotary wing aerodynamics in which t re tre a finite number of wakes below the rotor 
blades. It will be shown that this tin ry s completely analogous to Theodorsen and 
Loewy's work, and applicable to Gar ck': work. 

B. SCOPE 

Chapter II is a short review of tin wc k of Theodorsen [Ref. 1] and Garrick [Ref. 
3]. The material is presented in a rn me to allow one to understand the key 
assumptions made in each theory u about being overcome by the details. Chapter El 
contains an analysis of Loewy's work [Ri f. 2] and applies it to Garrick. The closed 
form finite wake theory is also presented here, as well as the special case of a single 
wake. Chapter IV is a limited discussion on applications of finite wake rotary wing 
unsteady aerodynamic theory. The fust section discusses logical extensions of this 
theory to include forward airspeed [Ref. 4] and compressibility effects [Ref. 5 and 6], 
A comparison of this theory to the Naval Postgraduate School (NPS) Unsteady Panel 
Code [Ref. 7 and 8] for the case of pure plunge is also included in the first section. 
The second section presents an explanation for the performance benefits due to Higher 
Harmonic Control (HHC). Appendix A is another application of this theory describing 
the efficient nature of flight of the hummingbird. The discussions in Chapter IV and 
Appendix A are not designed to provide all the answers to the questions arising from 
these applications. They are designed tc show that positive performance effects can be 
found in the application of rotary wing i nsteady aerodynamics. 
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II. BACKGROUND 



A. THEODORSEN'S LIFT DEFICIENCY FUNCTION 

The aerodynamic forces and moments on an oscillating airfoil for the fixed-wing 
case as determined by Theodorsen [Ref. 1] are based on potential flow theory and the 
Kutta condition. Potential flow theory reduces to Laplace's equation for the case of 
incompressible subsonic flow. 



d 2 $ + c?± x Q 



( 1 ) 



8x 2 dy 2 

Since potential flow theory follows the principle of superposition, the potentials can be 
treated as two types: non-circulatory and circulatory. The non-circulatory terms are 
primarily from the effects of the freestream with the airfoil. The circulatory terms 
result from the vorticity generated in the wake of the airfoil. It is in the circulatory 
terms that the lift deficiency function is defined. 

The aerodynamic forces and moments on a thin airfoil, depicted in Figure 1 , are 
obtained by integrating the pressure difference across the airfoil determined from the 
generalized Bernoulli equation as follows: 



P = 






b j p{x)dx 



( 2 ) 
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All distances are non-dimensionalized by the semi-chord (b) 



Figure 1 . Parameters of the thin airfoil. 



and 



AC ■ t > 2 J ‘ 



p(x)(x-a)dx 






where the generalized Bernoulli equation is written as 



^ - Mi* * f ) ■ 



( 4 ) 



In the manner suggested by Scanlan and Rosenbaum [Ref. 9], the potential flow 
for the non-circulatory terms can be written as 
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( 5 ) 






The pressure difference across the airfoil for the non-circulatory terms is 



A P(x)jv = 



-2pK r/n 



TT\f T-X 
/top b 

77 



/; 
H 






1 -;rl; - y^l — jr 2 \/ 1 -£ 2 

- + v/T-a^/T-^ 



( 6 ) 






where the downwash 1 w(£) = /j + b(£-a)<x + va, and simple harmonic motion is 
assumed, or w = w 0 e'“'. Integrating the pressure difference in equation (6) to obtain 
the non-circulatory forces and moments yields 

P N = -Trpb 2 (va + h- baa) (7) 



and 



(A/ = -7rp£ 2 [(— -a) vbo. + £ 2 (— + ^ 2 )a . 

1 2 8 



( 8 ) 



The potentials of the circulatory terms are determined by integrating the velocity 
potential resulting from a single vortex element with a vortex distribution per unit 
length, U(xJ, that is assumed to vary sinusoidally. This yields a circulation potential 
of 



‘Only two degrees of freedom (pitch and plunge) are presented here since most rotary 
wing aircraft do not have ailerons incorporated on the rotor blades. 
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( 9 ) 



r% oo /» i 

<P r = J </(p r = J 



U{x c )b 

-tan 1 



2xr 






1 - 






The Kutta condition, which requires the velocity at the trailing edge to be finite, can 
be written as 



(^) J^E) = finite. 

' dx 'edge ' dx 'edge 



( 10 ) 



This results in 



w(l) 



\ 



. -if' 

1-5 2 J, 



Jr + 1 

0U,)-* dx. . 



fiFi 



(ii) 



Now defining a lift deficiency function as 



l 



l <? -/A' A’ , 

e 



a*) = 






f 



*o + l -ikx a j 

e °dx„ 



\W~i 



02 ) 



the Kutta condition becomes 






JUTcu)- 



dx„ 



03) 



Using equation (4) the pressure difference across the airfoil for the circulatory tern* is 
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*p(x\ = 






ii£ 

-,N 1-5 






(14) 



Integrating the pressure difference in equation (14) to obtain the circulatory forces and 
moments yields 

P T = -2TrpvbC(k)\va + .A + b(^ -a)or] (15) 

and 

(A f) r = 277-p vb 2 (a+ C{k)\va + A+ . (16) 

The total pressure difference is obtained by adding the pressure differences for the 
non-circulatory and circulatory terms. After some algebraic manipulation the total 
unsteady pressure distribution can be written as 



a/? Or) 

pv 






l -a- pTT 

i+*N 1-5 






1 -X 



1+5 



77 



IN i \-i\x-i. 









(17) 



»K5)</5 



Likewise the total forces and moments are obtained by adding their two parts together 
yielding 



P - -Trpb 2 [ya + h- baa} - 2vp vbC{k)\^va + h+ />(— - aja 



(18) 
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and 



M a = -7rp/> 2 ^ -a) vb‘ + L ’(-i + <? 2 )a -abh]^ 

+ 277 p vtr d 4 - + .A + -a)aj . 



(19) 



The analytical solution to the lift efi iency function involves the use of complex 
Bessel functions (or Hankel functions , a; d can be written as 



C{k) 



H, 



( 2 ) 



H 2) ♦ iH™ * 



( 20 ) 



where Hj 2> = J. - iY n is the Hankel funct' >n of the second kind of order n. The lift 
deficiency function can be separated into its real and imaginary parts as 

C(k) - F(k) h iG(k) , (2D 



where 

m . W’*;)* 7,(y r / 0 ) 
Ur 

and 



( 22 ) 



G(k) = - 






(23) 



( y 1+ r 0 ) 2 + ( r ,-/ 0 ) 2 

and all Bessel functions are evaluated at the reduced frequency ( k ). A typical plot of 
F(k) and G(k) as given by Theodorsen [Ref. 1], is shown in Figure 2. Figure 3 is a 
semi-logarithmic plot of Theodorsens lift deficiency function, and the limits of F(k) 
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and G(k) as k approaches zero and infinity can be more easily seen. It is an 

interesting note that when = 0 , k = G = 0.188773655.... 

d k 

B. GARRICK S PROPULSIVE FORCE EQUATION 

One of the earliest investigations on the effect of an oscillating airfoil on the 
resultant drag was conducted by Katzma? r [Ref. 10] in two series of experiments. In 
his second series, Katzmayr measured a i egative drag or propulsive force for certain 
velocities and frequencies. This discovery, later to be known as the "Katzmayr 
effect", helped answer questions or how birds propel themselves through the air. 
Richardson [Ref. 11] used the "Katzmayr effect" to explain the locomotion of fish 
through the water. Garrick [Ref. 3] applied the equations of motion determined by 
Theodorsen [Ref. 1], equations (18) and (19), to a method developed by von Karman 
and Burgers [Ref. 12] to obtain a closed form analytical solution for the propulsive 
force generated by an oscillating airfoil. 2 

Garrick determines the average horizontal force by two methods: energy and 
force. He then compares the results of these methods, and shows they are identical. 
For simplicity only the force method is presented here. First he assumes simple 
harmonic motion for each of the degrees of freedom, resulting in 

a = a 0 e i{wt '^ (24) 



2 Though Garrick developed his equations for three degrees of freedom, only two 
degrees of freedom (pitch and plunge) will be presented here. 



10 



and 



h = h 0 e Ha, * 9l) , ( 25 > 

where (p 0 and (p 2 are phase angles between the initiation of the input and an arbitrary 
reference point. 

The total horizontal force on the airfoil is given by 

P x = npS 2 +aP (26) 



where P x is the propulsive force on the airfoil (positive to the left), P is the resultant 
force on the airfoil (positive down), and S is the component of induced downwash in 
the x-direction. The time averaged propulsive force is found by 

- (27) 

' 2tt Jo 

The x-component of the induced downwash is 



S 



-p|2 C(Xr)jva + A + b[^ - - Z)d| . 



(28) 



In finding the solution to the time averaged propulsive force equation, it is 
arbitrary whether one chooses to employ the real or imaginary parts; the results will be 
the same. Garrick [Ref. 3] chose the imaginary parts, and determined the time 
averaged propulsive force to be 
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( 29 ) 



P x = irp^ooM h 2 a [F 2 + G 2 ) 



b 2 a 



A [(-^(^ 2 * G 2 ) + J]t + -j) Sin((p 2 - <Po) 

*((| * ° 2 y f -^)cos«p 2 - (Po)| . 



It is interesting to note that in the special cases of only one degree of freedom, the 
effects of the arbitrary phase angles between initiation and the reference point 
disappear. The propulsive force equations become for pure plunge (a 0 = 0) 

P x = Trpbu 2 hl(F 2 + G 2 ) , ( 3 °) 



and for pure pitch (h 0 = 0) 3 



= npbu 2 b 2 a 2 0 \(F 2 + tf 2 )[-L + (j - a ) 2 ] + ± (j - a) 



k 2 '2 



a + 



2 '2 
LUl + a \ G 



k 2 ' ^2 +a > k 



n) 



( 31 ) 



If it is assumed that pitch and plunge are in phase with each other (they reach their 
maximum deflections at the same time), the coupled propulsive force equation 
becomes 



3 The result presented here differs from Garrick's equation (34) in reference 3. When 
equation (29) in reference 3 is simplified using the definitions of a 2 and b 2) the result is 
equation (31) above. This result has been confirmed by Garrick's work in reference 10. 
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( 32 ) 



P x = TTP^CD 2 + 

+ + ^ 2 )(-l- + (A _> 2 ) + A (A _) _ _ ^ _ (A ^ ] 



ba o h o 






Since propulsive force is analogous to negative drag, it is convenient to define a 
propulsive force coefficient as 



C p. ~ 



\pv 2 (lb) 



p v 



2 b 



( 33 ) 



The propulsive force coefficients for the special cases become for pure plunge 

(«. = 0 ) 

C Pm = TTk 2 h 2 0 (F 2 + G 2 ) , ( 34 ) 



where h 0 is the amplitude of deflection non-dimensionalized by the semi-chord, for 
pure pitch (h Q =0) 




and for coupled pitch and plunge 






( 36 ) 
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Ill / NALYSIS 



A. LOEWY S LIFT DEFICIENCY FL ECTION 

The aerodynamic forces and mor nt 1 - on an oscillating airfoil for the rotary-wing 

case are more complex than the fixed .vii \ counterpart. Loewy [Ref. 2] determines 

these forces and moments for the ho\ r c se where k = , but uses a slightly 

ftr 

different approach to the problem by solv. ig the integral downwash equation which 
leads to an equation for the pressure distr; xition which is in the same form as 
Theodorsen. Loewy then defines a modif ed lift deficiency function, and states that 
the integration of the pressure distribution across the airfoil would be identical to 
Theodorsen. This makes the expressions for forces and moments identical except for a 
modified lift deficiency function. 

The first step that Loewy performs i s to set up a system to account for the wakes 
generated by the previous blades in the same revolution and all blades in previous 
revolutions as shown in Figure 4. He uses two indices to account for the vorticity 
shed by a given wake: n indicates the re\ olution and q indicates the blade whose 
wake it is. The induced velocity or downwash resulting from an element of vorticity 
is obtained from the Biot-Savart Theorem, 



dw(x) = 



Y 



2 rr[(A - O 2 + (/?<?+ <7) 2 ^ 2 ] 



(37) 
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where y*, is the vorticity, Q is the number of rotor blades, and h is the non- 
dimensional wake spacing. Writing the integrals involving the bound vorticity and the 
vorticity in the wake of the reference airfoil separately from the rows of vorticity 
below the plane of the rotor yield 

i Ir' uodz , r- Yoo (g)</g , ^ " r v 

277 U-, x-l J, x-l h h> J.oc {x-D 2 + (nQ+q) 2 h 2 (3g) 

A r yJKHt-Ddl ' 

D= 1 j- oe (a~ 0 2 + -^ 2 Q 1 h 1 

In essence the problem has been broken down in the same manner as Theodorsen. 
The first integral represents the effects of the freestream on the airfoil (non-circulatory 
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terms). The second integral represents the velocity created by the vorticity generated 
by the reference wake (circulatory term). The third and fourth integrals represent the 
velocity created by the vorticity generated by previous blades or in previous 
revolutions (circulatory terms). The main difference between Loewy and Theodorsen 
is the terms which account for the vortici y generated by previous blades or in 
previous revolutions. 

Loewy [Ref. 2] shows that the vortici y shed by the qth blade in the mh revolution 
is 

Y = ikT e L '^ tq i t- 2nm 9lQ’ 27rmo '>] t (39) 



where T is the total circulation around the airfoil, vj/ 9 is the phase angle by which the 
motion of the <?th blade leads that of the reference blade, and m is the ratio of 
oscillatory frequency to rotational frequency. Substituting the vorticity expression 
from equation (39) into the integral downwash equation (38) yields 



•« • -± 






1 y,«¥5 






-m 



dl 






Q-\ 



e -i(2naqlQ-<i> g ) ^ ^ -rtirwn 
q = 1 n = 0 



r 



e ,k ^{x- £)</£ 

(x- i;) 2 + {nQ 4 of ti 1 



(40) 



IJcTYl e~ i24Tnw 



D ~ 1 



f 



(x - £) 2 + n 2 Q 2 h 2 - 



The last two integrals in equation (40) have the form 



e'^'jx-Jj') a 'j _ J - v . e -mx*A) 
« {x-tf + A 2 



( 41 ) 
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Substituting equation (41) into equation (40) and noting that the summations over n 
are convergent geometric series give 




where 



W{k,hjri) = 



^2 (^ e ^Q e i2nwjQ-</)l(? e ^ 9 
Szl 



( 43 ) 



gkhQ gilTrm _ j 



The function W(k,h,m) may be thought of as a weighting function for the vorticity 
shed by previous blades or in previous revolutions. 

The form of the downwash equation in equation (42) can be solved by applying 
Sohngen's inversion formula [Ref. 14], which shows that the solution to an equation in 
the form 



Satisfying the condition f(l ) = finite is the same as employing the Kutta condition. 
The bound vorticity becomes 




( 44 ) 



is 
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Evaluating the circulation over the entire airfoil yields 



T = e ik 
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Y Ax) dx = 
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where the Hankel and Bessel functions are evaluated at reduced frequency ( k ). 

Since the airfoil can be thought of as a vortex sheet, the generalized Bernoulli 
equation (4) becomes 



a p(x) = -p 






( 48 ) 



If simple harmonic motion is assumed, 
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and equation (48) becomes 



Lp(x) = -p 



yy a (x,t) + /o )bj Y a (U)dl 



( 50 ) 



Substituting the bound vorticity equation (46) and the airfoil circulation equation (47) 
into equation (50) yield 
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The pressure distribution in equation (51) has the same form as Theodorsen's in 
equation (17) if the factor multiplying the first integral is written as 




C\kji,m)\ , 



(52) 



where C'(k,h,m ) is Loewy's lift deficiency function. Solving for C'(k,h,m) yields 



CXk.hyin) = 
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(53) 



where the Hankel and Bessel functions are evaluated at reduced frequency (k). Since 
C'(k,h,m) is not a function of location along the airfoil, the integration of the pressure 
distribution in equation (51) across the airfoil will yield equations of motion identical 
to Theodorsen's in equations (18) and (19) except C(k) will be replaced by C'(k,h,m). 

It can be seen that as W(k,h,m) approaches zero, C'(k,h,m) = *7 (it), and this condition 
corresponds to an infinite wake spacing (/z-~). 

Since the wake weighting function is periodic in nature, Loewy [Ref. 2] shows that 
the wake weighting function for a multi-blade rotor can be expressed by a single-blade 
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rotor with modified values of h and n tli yield the same value of W. For a single- 
blade rotor, the wake weighting func' >n >ecomes 



WikJijn' 
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(54) 



where h and m are the modified vah.. .. o: wake spacing and frequency ratio 
respectively. For the case of the sir g c-fc ide rotor the wake weighting function can 
be separated into its real and imagin / p rts, 



W(k L-n) = a + $ , 



(55) 



where 



and 



cos 2 rrm - e^_ 

e kh Icoslvm + e~ kb 



(56) 



$ = * in 2nm . (57) 

e kb - 2 cos 2tt m + e~ kb 

Separating Loewy's lift deficiency function into its real and imaginary parts similar to 
Theodorsen yields 4 

C'ikjKw) = F' + iG’ ( 58 ) 



where 



4 This is equivalent to Loewy's equations in reference 2 when the substitutions 
a = 1 + 2a and P = -2(3 are used. 
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( 59 ) 
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and 
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Plots of Loewy's lift deficiency function are shown in Figures 5 through 8. It is easily 



seen that there are rapid changes in the lift deficiency function as reduced frequency 
(k), wake spacing (h), and frequency (m) vary. 

Since Loewy's lift deficiency function is completely analogous to Theodorsen's lift 
deficiency function, the propulsive force determined by Garrick [Ref. 3] and the 
propulsive force coefficient can be calculated for the rotary-wing case by substituting 
F' and G' for F and G. For pure plunge (a 0 = 0) 




(63) 



for pure pitch (h 0 = 0) 



C p = 7T 
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(64) 
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Figure 5. Loewy's lift deficiency function as a function of wake spacing for 

m = 0 . 
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Figure 6. Loewy's lift deficiency function as a function of wake spacing for 
m = 0 . 25 . 
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Figure 7. Loewy's lift deficiency function as a funciton of wake spacing for 
m = 0.5. 
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Figure 8. Loewy's lift dieficiency function as a function of wake spacing for 
m - 0 . 75 . 
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and for coupled pitch and plunge 
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Plots of propulsive force coefficient a a 5 inction of frequency ratio for = 0.1234 are 
shown in Figure 9 through 11. It c< 1 be : een that Cpx is reasonably constant for 0.15 
< m < 0.8 and greater that the "Katzmnyr effect." (In this case h = 100 is sufficient to 
represent /?-«>). For m < 0.15 and m > 0.8, Cpx begins to decrease to the point where 
it eventually becomes less than the "Katzmayr effect." 

B. FINITE WAKE LIFT DEFICIENCY FUNCTION 
1. General 

Loewy's theory can be extended to the case of a finite number of wakes 
below the rotor blade by modifying the wake weighting function. This modified wake 
weighting function is completely analogous with Loewy's which makes the definition 
of the lift deficiency function the same, but with a and $ modified by the new wake 
weighting function. Since the lift deficies cy function is of the same form, the 
aerodynamic forces and moments will be the same, but with a modified lift deficiency 
function. Loewy has already shown that he rotary- wing problem reduces to the case 
of a single rotor. The modified wake w eighting function will also be developed for 
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Figure 9. Propulsive force coefficient in pure plunge as a function of wake spacing 
with an infinite number of wakes (Loewy's lift deficiency function). 




Figure 10. Propulsive force coefficient in pure pitch as a function of wake spacing 
with an infinite number of wakes (Loewy's lift deficiency function). 
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Figure 11. Propulsive force coefficient in coupled pitch-plunge as a function of wake 
spacing with an infinite number of wakes (Loewy's lift deficiency function). 



the case of the multi-blade rotor, and then reduced to the case of the single-blade, rot or 
with modified wake spacing and frequency ratio. 

The development of a finite wake lift deficiency function is the same as 
Loewy's described earlier. Changing the integral downwash equation (40) to a finite 
number (AO of revolutions yields 
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where the summations from n equal zero to infinity have been replaced by finite 
summations from n equal zero to N. Using equation (41) to solve the last two 
integrals in equation (66) results in 



w{x) 



277 r 



-i *-5 




e ~ ik l d X 

x-l 



. Q-\ N 

jj-JcT e Vy; e -i2TTmn e ~k{nQ+q)b 

' ?= 1 D-0 

N 



+ ^ -i 'Zimin g -nQkh \ 

i7 — 1 ■ 



(67) 



Since this is in the same form as equation (42), the terms enclosed by the parentheses 
can be defined as the finite wake weighting function, or 



Wj^k,hyin) = 
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It is easily seen that this modified wake weighting function is periodic once the 
relationship i|; (? is known, and a multi-blade rotor can be reduced to a single-blade 
rotor with modified wake spacing and frequency ratio. For the case of the single- 
blade rotor {Q = 1), 



WjJUcJi,w) = 
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g-ilrrmn g-nkh 
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(69) 



It can be shown that as N-oo, the finite wake weighting function in equation (69) 
becomes a convergent geometric series equivalent to Loewy's wake weighting function 
in equation (54) provided m * 0 and h * 0. A singularity in the solution exists when 
m = 0 and h = 0, and the case of h = 0 has no real physical significance. 
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The solution to the remainder of the finite wake problem is identical to 
Loewy, but using the modified finite wake lift deficiency function. To keep the finite 
wake lift deficiency function separate from Loewy's lift deficiency function, the former 
is defined as follows: 
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Separating the finite wake weighting function in equation (69) into real and imaginary 
parts, 
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yields 
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The finite wake lift deficiency function becomes 

C' = F'+iG’ , 



(72) 
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where 
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(75) 
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and 
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2. Single Wake Case 

For the case of a single wake, the wake weighting function reduces to 

W | = e -(i2nm + kh) 



(79) 



where 



a, = e ~ kb cos 2Trm 



(80) 



and 
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Plots of the single wake lift deficient:; f*. ction are shown in Figures 12 through 15. 
The propulsive force coeffic :it an be calculated by replacing F and G in 



equations (34), (35), and (36) with / tnc G\ For pure plunge (a„ = 0) 



c Pj - nr- ;{ : f-) 2 +(G') 2 ) , 



( 82 ) 



for pure pitch (h 0 = 0) 
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and for coupled pitch and plunge 



C P = 7rk 2 \h 2 0 UF-) 2 + (G') 2 \ 



2 '2 1 h k 2] '2 > k 

«x\(\ - *) (ir )• ■• (c- yy \ 4 - - (f- ) ) 



( 84 ) 



Plots of the propulsive force coefficient as a function of frequency ratio for k = 0.1234 
are shown in Figures 16 through 18. It is interesting to note that for 0.24 < m < 0.7 
the propulsive force is greater than the "Katzmayr effect" with the maximum occurring 
just prior to m = 0.5. Comparing Figure 9, 10, and 11 with Figures 16, 17, and 18 
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Figure 12. Single wake lift deficiency function as a function of wake spacing for 

m = 0 . 
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Figure 13. Single wake lift deficiency function as a function of wake spacing for 
m = 0.25. 




Figure 14. Single wake lift deficiency function as a function of wake spacing for 
m = 0.5. 
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Reduced Frequency (k) 



Figure 15. Single wake lift deficiency function as a function of wake spacing for 
m - 0.75. 
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Figure 16. Propulsive force coefficient in pure plunge as a function of wake spacing 
with a single wake. 
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Figure 17. Propulsive force coefficient in pure pitch as a function of wake spacing 
with a single wake. 
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Figure 18. Propulsive force coefficient ii coupled pitch-plunge as a function of wake 
spacing with a single wake. 

respectively, the maximum propulsive force (at approximately m = 0.5) for a single 
wake is always greater than the infinite number of wakes. The reason for this increase 
is that the absolute value of C and C re never greater than 1.0, but C* as m 
approaches 0.5 exceeds 1.0 for certain wake spacings. 

When m = 0.5, this corresponds to the wakes being 180° out of phase as 
shown in Figure 19, and the vortex from the wake of the reference blade in the current 
revolution is directly above an oppositely spinning vortex from the single wake of the 
previous blade or revolution. This efficient use of the vortices in the wake causes C * 
to exceed 1.0. As the number of wakes increases the vortices from lower wakes are 
not always 180° out of phase with the reference blade. When m - 0.5, each wake is 
180° out of phase with the wake above a td below it. Thus, each wake is in phase 
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with the wakes ±2 wake spacings from it. Since the finite wake weighting function 
contains the term e' ,2nm , the value of each term in the wake weighting function will 
alternate in sign. However, since the wake spacing is greater for each subsequent 
wake, the absolute value of each term is not as large as the first term, and thus, the 
effect of the in-phase wakes are not as strong. As the wake spacing increases this 
effect diminishes. 

3. Multiple Wake Case 

When the number of wakes increases, the value of frequency ratio at which 
the maximum propulsive force coefficient occurs changes as shown in Figures 20 
through 22 for h = 2.0. This change is analogous to the estimation of a square wave 
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by a finite number of terms in a Fourier series, where the propulsive force coefficient 
determined using Loewy's wake weighting function is considered the square wave. As 
Figures 16 through 18 and 20 through 22 show, an increase in either wake spacing ( h ) 
or the number of wakes (AO will decrease the maximum propulsive force that can be 
obtained. The limit of this decrease will be Loewy's case for a hover. 




Figure 20. Propulsive force coefficient in pure plunge as a function of the number of 
wakes at h = 2.0. 
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Figure 21. Propulsive force coefficient in pure pitch as a function of the number of 
wakes at h = 2.0. 




Figure 22. Propulsive force coefficient in coupled pitch-plunge as a function of the 
number of wakes at h = 2.0. 
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IV. APPLICATIONS 



A. EXTENSIONS AND COMPARISON OF THEORY 
1. Forward Airspeed 

Extension of the finite wake theory can easily be made to include forward 
airspeed. Shipman and Wood [Ref. 4] expand on Loewy's work incorporating forward 
airspeed via the advance ratio (p) into the rotary wing unsteady aerodynamic problem. 
The definition of reduced frequency becomes 



k = 



(1 ‘ 



(85) 



The Biot-Savart Theorem is used with wakes extending out to infinity. Since 
freestream velocity at any point on the rotor blade is a function of both blade section 
radius and azimuth position, the shed vorticity will experience a build-up and decay as 
azimuth varies from 0° to 360° . Shipman and Wood developed a decay function to 
account for the variation in vorticity and incorporate this function and the advance 
ratio into the integral downwash equation. The pressure distribution is found by 
applying Sohngen's inversion formula [Ref. 14] to the downwash equation. The 
pressure distribution is shown to be in the same form as Theodorsen and Loewy, and a 
modified lift deficiency function is defined to account for the differences. Since this 
modified lift deficiency function includes Loewy's wake weighting function as one of 
its terms, the finite wake weighting function could be directly substituted into Shipman 
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and Wood's lift deficiency function. If the segment of the wake weighting function 
due to the vorticity decay (AW) is modified to account for a finite number of wakes 
instead of an infinite number, a new finite wake forward airspeed lift deficiency 
function can be defined, and the remainder of the problem is the same as Shipman and 
Wood. 

2. Compressibility Effects 

Adding the effects of compressibility to the rotary wing unsteady aerodynamic 
problem was first accomplished by Jones and Rao [Ref. 5]. Hammond [Ref. 6] 
provided an alternate approach to the problem shortly after Jones and Rao, but only 
the method by Jones and Rao will be presented here. 

Jones and Rao expand on Loewy's work by taking the known solution to the 
fixed wing two-dimensional unsteady airfoil theory and modifying it to include the 
effect of an infinite number of wakes below the reference airfoil (or rotor disk). 

Euler's equation is used with a coordinate change to determine the downwash equation. 
An interesting result is found, and that is for a given circulation, the downwash due to 
an infinite system of wakes is the same for compressible and incompressible flow. In 
other words, the circulation generated by layers of shed vorticity beneath an airfoil is 
independent of Mach number. This fact will be true whether there is an infinite 
number of wakes or a finite number of wakes. Thus, the finite wake theory can be 
applied here, and the remainder of the problem is the same as Jones and Rao. 
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3. NPS Unsteady Panel Code 

The NPS Unsteady Panel C~ ; v as written by Pang [Ref. 7] to solve the 
potential flow for two airfoils execut . j u isteady motion in an inviscid, 
incompressible flow. In this code, tl. tw - airfoil surfaces are approximated by a large 
number of surface elements, and a un or i source distribution and vorticity are placed 
on each element. The source streng 1 vr es from element to element , while the 
vortex strength is the same for all e.e - en ;. The singularity strengths are determined 
from the flow tangency condition oi. ot) airfoil surfaces and the two Kutta 
conditions. 

Extensive comparisons of th.- co le with Theodorsen was accomplished by 
Riester [Ref. 17], and the code (using an NACA 0007 airfoil) shows very good 
agreement with flat plate theory. Compa'ison of this code with Loewy could not be 
accomplished easily since it would require a very large number of airfoils to generate 
the infinite number of wakes. Using ti e finite wake theory, the code could easily be 
compared to the case of a single wain . A comparison for pure plunge [Ref. 8] with h 
= 2.0, k = 0.1234, and h 0 = 0.14 are presented in Figure 23. The results show that the 
NPS Panel Code has very good agreement with the finite wake theory. 

B. HIGHER HARMONIC CONTROJ 

Higher Harmonic Control (HHC) is an active vibration control system for 
helicopters. The concept of HHC is to ader the aerodynamic loads on the rotor blades 
such that the blade response is reduced, vhich in turn reduces the vibratory forces and 
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Figure 23. Comparison of propulsive force coefficient in pure plunge (h = 2.0, k = 
0.1234, and h 0 = 0.14). 

moments acting at the hub that cause vibration. HHC flight test results [Ref. 15], 
under a joint NASA-Army sponsored program using a modified OH-6A, show a 
reduction of airframe vibrations up to 90 percent. Basically, HHC is an electronic, 
computer-controlled active vibration suppression system which senses and cancels 
vibrations in a helicopter by Q per revolution feathering or pitch motion of the rotor 
blades, where Q is the number of rotor blades. 5 On the OH-6A, higher harmonic 
blade pitch control was achieved by superimposing 4/rev (or 4P) swashplate motion on 
basic cyclic and collective control inputs. The aircraft was flown from hover to 100 
knots, but the remainder of the discussion will focus on the hover performance. 



5 Normally for helicopters, N is used to represent the number of blades, but to 
remain consistent with Chapter III Q will be used. 
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In addition to reducing vibration levels HHC on the OH-6A showed the potential 
for decreased helicopter power requirements. Figures 24 and 25 show hover power 
required for three different schedules of HHC open loop excitation. What is of 
interest here is that the amount of power reduction achieved (on the order of 10% or 
20 horsepower) appears to be relatively independent of the type of excitation applied. 
This power reduction is important, especially in hover where power required is large, 
because of the benefit of lowering fuel consumption, which can be translated into a 
larger payload or increased mission time. 

The physics behind HHC can be explained by unsteady aerodynamics. The input 
by HHC to the rotor blade is a pitch oscillation, ±Vz degree on the OH-6A. From 
Figure 10 (infinite wakes) it can be seen that the propulsive force coefficient is 
negative (drag) for the wake spacing of the OH-6A (h - 2.0) regardless of the phasing 
of the wakes (m). Thus, it might be concluded that pitch oscillations alone cannot 
account for the substantial reduction of power recorded by the flight tests. If it is 
assumed that viscous interaction of vortices reduce the problem to a finite wake 
solution as shown in Figures 17 and 21, the propulsive force obtained is still not 
sufficient to account for the power reduction even if the wakes were optimally phased. 

As stated earlier, HHC provided a 4P pitch input in both collective and cyclic were 
applied at the blade root. The collective input results in a 4P oscillation along the 
blade, while the cyclic input results in 3P and 5P oscillations along the blade. From 
the Southwell plot of the OH-6A [Ref. 16] shown in Figure 26, it can be seen that the 
second flapwise and third flapwise natural frequencies occur near 3P and 5P 
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SHAFT HORSEPOWER 




Figure 24. Effect of HHC on main rotor shaft power in hover for selected set of three 
open loop inputs. 

respectively. During the course of a revolution of a blade about the hub, the blade 
will naturally experience flapping motion due to changes in angle of attack, or pitch 
changes, as the blade proceeds around. Since 3P and 5P flapping frequencies are 
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SHAFT HORSEPOWER 



AIRSPEEO-HOVER 




Figure 25. Effect of HHC on main rotor shaft power in hover for selected set of two 
open loop inputs. 

near natural frequencies, small 3P and 5P pitch inputs will result in much greater 
flapping outputs. Thus, the unsteady aerodynamics of HHC for the OH-6A is coupled 
pitch-plunge oscillations. Figures 11, 18, and 22 show that coupled pitch-plunge 
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FREQUENCY RATIO (CYC/REV) 




Figure 26. Southwell plot of the OH-6A rotor blade. 
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V. CONCLUSION 



The unsteady aerodynamics of rotary wings is inherently more complicated 
because of the layers of shed vorticity, or wakes, beneath the rotor disk that interact 
with the rotor blades. The Biot-Savart Theorem is used to account for the effects of 
these wakes, and a modified lift deficiency function is developed for use with 
Theodorsen's equations of motion [Ref. 1]. Loewy [Ref. 2] develops a lift deficiency 
function for the case of an infinite number of wakes, while this thesis develops a lift 
deficiency function for any finite number of wakes. This thesis also presents the 
special case of a single wake, and shows that for small wake spacings there are certain 
reduced frequencies where the lift is actually enhanced. In other words, the lift 
deficiency function becomes a lift efficiency function. 

The equations describing the propulsive force generated by an oscillating airfoil 
are developed by Garrick [Ref. 3]. These equations are derived from Theodorsen's 
equations of motion, and included the lift deficiency function as a parameter. Since 
the rotary wing lift deficiency functions (Loewy's and the finite wake) are just 
modifications of Theodorsen's work, Garrick's work can be applied to rotary wing 
unsteady aerodynamics. When either Loewy's lift deficiency function or the finite 
wake lift deficiency function is used, the propulsive force can be greatly enhanced 
with the proper phase relationship of the wakes. The use of coupled pitch and plunge 
can result in a propulsive force which is an order of magnitude higher than either pure 
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plunge or pure pitch. This coupling of plunge and pitch is used to explain the 
performance benefits reported in the flight test data for the OH-6A equipped with 
Higher Harmonic Control. Coupled pitch-plunge oscillations with proper phase 
relationships of the wakes yield a propulsive force large enough to account for a 
substantial power reduction for the helicopter. 
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APPENDIX A 



I. FLIGHT OF THE HUMMINGBIRD 

For thousands of years, man has been intrigued by the flight of birds. Early Greek 
mythology includes a tale of Daedalus and Icarus who made wings held together by 
wax, and escaped their imprisonment on the island of Crete by flying through the air 
like birds. Early attempts by man to fly were usually designed around imitations of 
birds. Most of these early attempts met disastrous results because of a lack of clear 
understanding of the aerodynamics of ornithological flight, and that is the wing on a 
bird is a coupled lift and propulsion device. In fact, the Wright brothers first 
successful flight demonstrated man's limited knowledge of aerodynamics by 
decoupling lift and propulsion in order to achieve flight. To this day lift and 
propulsion remain decoupled. 

The coupling of lift and propulsion in birds can be explained by unsteady 
aerodynamics. As a bird flaps its wing down, the trailing edge also begins to rise 
[Ref. 18]. Thus, the motion of the wing is a coupled pitch-plunge motion. It was 
shown earlier that coupled pitch-plunge motion yields a propulsive force that is an 
order of magnitude higher than pure plunge. Since most birds require some forward 
airspeed before being capable of achieving flight and the rate of flapping is between 1 
to 5 beats per second, the analogy to unsteady aerodynamics can be described by 
Garrick's equation for coupled pitch-plunge using Theodorsen' lift deficiency function 
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(equation (36)). Birds take advantage of the extra propulsive benefit of coupled pitch- 
plunge to make flight appear effortlessly. 

A few birds, most notably the hummingbirds, are capable of hovering. In hovering 
flight the hummingbird's body is angled upward at about 45°, so that the plane of 
wing flapping is approximately horizontal. As the wing moves up (backward) and 
down (forward), the angle of attack, or pitch, changes making a figure eight pattern 
that is on its side. In essence, the hummingbirds wing motion is similar to the motion 
of the hands of a swimmer when he is treading water. Once again the hummingbird 
takes advantage of a coupled pitch-plunge motion to achieve a higher propulsive force. 

What is unique about the hummingbird is the rate of flapping, or the oscillation 
frequency; it is on the order of 40 to 80 beats per second. Welty [Ref. 18] states that 
the hovering flight of the hummingbird is grossly inefficient when compared to 
helicopter rotor blades, but this is overcome by an increase in apparent velocity 
created by trailing currents from the previous reverse stroke. This is not exactly 
correct. The hummingbird's wing motion is similar to that of helicopter rotor blades. 
Both are cyclic motions of changes in pitch (angle of attack) and plunge (flapping). 
The hummingbird is just as efficient as a helicopter - both have a large power required 
in hover. The apparent increase in velocity can be explained by trailing vortices. 

Since the hummingbird's rate of flapping is large in a hover, the wake created by a 
previous stroke will not dissipate and be carried downstream. Each stroke of the wing 
always has layers of shed vorticity directly below it. If these wakes are optimally 
phased depending on the number of wakes, the benefit will be a propulsive force 
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greater than that of non-hovering birds. The hummingbird must use wakes efficiently, 
or it will not hover. Thus, the analogy of the hummingbird to unsteady aerodynamics 
can be described by Garrick's equation for coupled pitch-plunge using the finite wake 
lift deficiency function, or in the case of a pure hover Loewy's lift deficiency function. 
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